Multiple Andreev Reflections in a Carbon Nanotube Quantum Dot 
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We report resonant multiple Andreev reflections in a multiwall carbon nanotube quantum dot 
coupled to superconducting leads. The position and magnitude of the subharmonic gap structure 
is found to depend strongly on the level positions of the single-electron states which are adjusted 
with a gate electrode. We discuss a theoretical model of the device and compare the calculated dif- 
ferential conductance with the experimental data, (pdf including figures, see: www.unibas.ch/phys- 
meso/Research/Papers/2003/MAR-MWNT.pdf) 



PACS numbers: 74.78.Na,74.45.+c,73.63.Kv,73.21.La,73.23.Hk,73.63.Fg 
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The electronic transport properties of quantum dots 
coupled to metallic leads have been the object of exten- 
sive theoretical and experimental study . When weakly 
coupled to its leads, the low-temperature transport char- 
acteristics are dominated by size and charge quantization 
effects, parameterized by the single-electron level spacing 
AE and the single-electron charging energy Uc ■ When 
the coupling of the quantum dot to the source and drain 
electrodes is increased, higher-order tunneling processes 
such as the Kondo effect become important 0. New ef- 
fects are expected when the leads coupled to the quantum 
dot are superconductors. In that case electron transport 
at small bias voltages is mediated by multiple Andreev 
reflection (MAR) 0, 0. Unlike conventional S-N-S de- 
vices, however, the MAR structure is now expected 00 
to strongly depend on the level positions of the single- 
electron states of the quantum dot which can be tuned 
with a gate electrode. The influence of Coulomb and 
Kondo correlations have been addressed theoretically in 
Refs. 0. 

Because MAR is suppressed rapidly for low- 
transparency junctions, its observation requires a rela- 
tively strong coupling between the leads and quantum 
dot. Even more so as on-site Coulomb repulsion, which 
is common to weakly coupled dots, is expected to re- 
duce Andreev processes even further. A quantum dot 
very weakly coupled to superconducting leads has been 
studied experimentally by Ralph et al. 0. In this case 
the transport characteristics were indeed dominated by 
charging effects and MAR was completely suppressed. 

The coupling to the leads, expressed in the life-time 
broadening T of the levels of the quantum dot, should 
be compared to the superconducting gap energy A. 
Favourable for the observation of MAR in a quantum 
dot are coupling strengths of order T ~ A and a small 
charging energy Uc < A. Together with the restriction 
that T < AE (for any quantum dot) this leads to the 
approximate condition Uc < T < AE. For most quan- 
tum dots typically the opposite is true and AE <C Uc- 
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FIG. 1: (a) Greyscale representation of the differential con- 
ductance as a function of source-drain (V„d) and gate voltage 
(V g ) at T = 50 mK and B = 26 mT for a MWNT quantum 
dot. Here and in the following greyscale plots, the darker a 
region the higher the differentical conductance. The dashed 
lines outline the onset of first-order tunneling processes, (b) 
Linear-response conductance G as a function of V g . The ap- 
pearance of a single broad peak is due to the Kondo effect, 
(c) Differential conductance at two different values of V g . 



It has recently been shown however, that well-coupled 
multiwall carbon nanotube (MWNT) quantum dots can 
have favourable ratio's of AE/Uc ~ 2 and Uc can be as 
small as 0.4 meV, comparable to the energy gap 2A of a 
conventional superconductor like Al. 

Here we report on the experimental study of resonant 
MAR in a MWNT quantum dot. The superconducting 
leads to the MWNT consist of an Au/Al bilayer (45/135 
nm). Before investigating the system in the supercon- 
ducting state, the sample is first characterized in the nor- 
mal state by applying a small magnetic field. From these 
measurements relevant parameters such as AE, Uc, and 
r are obtained. We then discuss a theoretical model that 
describes the differential conductance of an individual 
level in a quantum dot coupled to superconducting elec- 
trodes. In the final part of the paper we compare the 
calculated differential conductance with the experimen- 
tal data. 

Figure ^ shows a greyscale representation of the dif- 
ferential conductance dI/dV s d versus source-drain (V s d) 
and gate voltage (Vg) at T = 50 mK when the contacts 
are driven normal by a small magnetic field. The dotted 
white lines outline the onset of first-order tunneling and 
appear when a discrete energy level of the quantum dot 
is at resonance with the electrochemical potential of one 
of the leads. ,[From these and other electron states mea- 
sured for this sample (about 20 in total), we obtain an 
average single-electron level spacing AE ~ 0.6 meV and 
a charging energy Uc ~ 0.4 meV. Since Uc = e 2 /C% this 
yields Cs = 400 aF for the total capacitance which is the 
sum of the gate capacitance (C g ) and the contact capac- 
itances C s (source) and Cd (drain). From the data of 
Fig. HJi we obtain C g /C s = 0.0036, and C s /C d = 0.45. 
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FIG. 2: (a) Greyscale representation of the calculated dif- 
ferential conductance as a function of V s d and level position 
(ip g oc Vg) for a single-electron level coupled symmetrically to 
superconducting electrodes. For clarity, the low-energy part 
eVsd < 2A/4 has been omitted. The dashed lines indicate 
resonance positions, as explained in the text, (b) Differen- 
tial conductance at ip g = 0. Note, that MAR peaks at 2A/n 
are suppressed for even values of n. (c) Schematics of a sin- 
gle electron state between superconducting source and drain 
electrodes. The situation shown corresponds to point C in 
panel (a). The spectral density is shown at the bottom. The 
Lorentzian level broadening in the normal state is replaced by 
a narrow central resonance accompanied by a series of satellite 
peaks. (d) Same as (c) for point D in panel (a). 

The lifetime broadening T is obtained from the width 
of the single-electron peaks at finite source-drain voltage 
(taking the background conductance into account) and is 
found to be T ~ 0.35 meV. 

The high-conductance ridge around V s d = mV in 
Fig. nji is a manifestation of the spin-1/2 Kondo effect 
occurring when the number of electrons on the dot is 
odd. As a result, the Coulomb valley in the conductance 
has disappeared in this region of V g and at 50 mK, which 
is the base temperature of our dilution refrigerator, only 
a single peak remains, see Fig. ^d. This is known as 
the unitary limit of conductance. The appearance of the 
Kondo effect is an indication that the coupling to the 
leads is relatively strong. We will not discuss the Kondo 
effect here, and instead refer to Ref. [Toj . 

When the magnetic field is switched off the leads be- 
come superconducting. To calculate the expected differ- 
ential conductance in the superconducting state of the 
leads we have used the non-equilibrium Green-function 
technique We model the quantum dot as a series 

of spin-degenerate resonant levels coupled to supercon- 
ducting electrodes, which are assumed to have a BCS 
spectral density. Note, that neither electron-electron in- 
teraction (Coulomb blockade) nor exchange correlations 
(Kondo effect) are accounted for in the model, which 
may, therefore, not explain all details of the actual mea- 
surements. However, the interplay between MAR and 
resonant scattering already leads to strongly nonlinear 
IV-characteristics and reproduces some of the key fea- 
tures of the data. The details of the model calcula- 
tion will be presented elsewhere |l2j . The main parame- 
ters entering the calculation are the two tunneling rates 
r ' s (d) = ^\t s (d) | 2 -^0j related to microscopic tunneling 
amplitudes t s (d) via the density of states Nq of the source 
(drain) contact, fn the model we account for the gate 
voltage by a shift of the level, which can be adjusted 
according to the experimentally observed Coulomb dia- 
monds, see Fig. ^ 

The discrete nature of the single-electron states is most 
pronounced when T is small. Therefore, before present- 
ing the model calculation that directly compares to the 
experimental data, we first discuss the transport charac- 
teristics of a single spin-degenerate level with a relatively 



FIG. 3: (a) Greyscale representation of the measured differ- 
ential conductance as a function of V 3 a and V g at T — 50 mK 
with the leads in the superconducting state. The gate volt- 
age range corresponds to the left part of Fig. la. The dashed 
white lines emphasize the position of the MAR peaks, (b-c) 
Differential conductance at the positions given in panel (a). 



weak and symmetric coupling to the superconducting 
leads along the lines of Refs. 0, The total tunnel- 
ing rate r = T s + Yd is set to A. Figure |2 shows the 
corresponding greyscale representation of the calculated 
differential conductance dI/dV s d versus ip g :— eV g C g jC^ 
and eV s d- The peak structure in dI/dV s d at V s d < 2A/e 
is the result of MAR. fn general, Andreev channels be- 
come available for transport at voltages V s d — 2A/ne, 
where n is an integer number. These positions are indi- 
cated by the horizontal dashed black lines in Fig.|2^. The 
appearance and magnitude of the MAR peaks, however, 
is strongly dependent on the position of the resonant level 
in the quantum dot with respect to the Fermi energy of 
the leads. Only those MAR trajectories that connect the 
resonant level to the leads' BCS spectral densities give 
a significant contribution to the current. Consider, for 
example, the position marked by C in Fig. which 
corresponds to the schematics of Fig. |2t and indicates 
the position (ip g , eV s d) = (0, 2A/3). The corresponding 
second-order Andreev trajectory connects the gap edges 
of the source and drain electrodes and includes the res- 
onant level which is situated exactly in between the re- 
spective Fermi energies. This results in the large peak in 
dI/dV s d seen in Fig. EJa. 

A similar peak is absent at (<p g ,eV s d) — (0, A), corre- 
sponding to point D in Fig. At this voltage, first- 
order Andreev reflection becomes possible. The corre- 
sponding trajectories (see Fig. [SJi), however, do not di- 
rectly connect the resonant level to the leads spectral 
densities, and therefore do not significantly contribute to 
the current. Only when the gate voltage is adjusted to 
align the level with the Fermi energy of one of the leads 
(indicated by the arrows in Fig.^t) a peak in dI/dV s d is 
observed, ft can be shown (for symmetric junctions) that 
the subharmonic gap structure at V g = is suppressed 
for all voltages V s d = 2A/ne with n = even 0, When 
V s d is increased beyond A/e, peaks are observed either 
when the level stays aligned with the electrochemical po- 
tential of the leads (red dashed lines: ip g = ±eV S( j/2) or 
when the level follows the gap edges as an initial or fi- 
nal state of an Andreev process (blue dash-dotted lines: 
<p g = ±(A - eV s d/2) or y g = ±(A - 3eV sd /2)). 

We now turn to the actual measurements of the differ- 
ential conductance when the leads are in the supercon- 
ducting state. Figure |21 shows a greyscale representation 
of the measured dI/dV s d versus V s d an d V g at B = OmT 
for the same single-electron state of Fig. ^ A number of 
differences between the normal state (Fig. ^) and super- 
conducting state (Fig. [3Jl can be observed. The horizon- 
tal high-conductance lines at V s d = ±0.2 mV in Fig. |3 
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FIG. 4: (a) Greyscale representation of the calculated differ- 
ential conductance as a function of V s d and level position (<p g ). 
The different adjustable parameters represent the experimen- 
tal situation, see text, (b-c) Differential conductance at the 
positions given in panel (a) . For the two dashed curves, V has 
been chosen approximately twice as large as the experimental 
value. 



for example, are attributed to the onset of quasi-particle 
tunneling when V s d — 2A/e. The subgap structure at 
V sd < 2A/e (i.e. below 0.2 mV) is attributed to MAR. 
As anticipated, the magnitude and the position (dashed 
white lines) of MAR peaks depend on Vj. To allow for 
comparison with theory, the adjustable parameters of the 
model are set to the values obtained from the measure- 
ment of Fig. ^ The most important parameter is the 
coupling between the electrodes and the dot which turns 
out to be r ~ 3.5 A. The voltage division between the 
two tunnel barriers separating the quantum dot from the 
leads is C s /Cd — 0.45. The individual tunneling rates are 
not exactly known but are not expected to show a strong 
asymmetry The neighboring single-electron states, 
separated by AE ~ 6.5 A, are included in the calcula- 
tion. The finite temperature of the experiment is taken 
into account and set to T = 0.1 A. 

The resulting calculated greyscale representation of the 
differential conductance is shown in Fig. 21 The over- 
all appearance clearly resembles the measured data of 
Fig. |21 For example, both the model and the measured 
data show a large peak in dI/dV s d around V s d = OmV 
when the electron state is at resonance with the electro- 
chemical potential of the leads (i.e. at V g — 0). When 
the level is moved away from this position, the linear- 
response conductance rapidly decays to values below its 
normal-state value. In contrast, the differential conduc- 
tance peak at V s d — 2A/e shows the opposite behavior 
(both in the model as in the experiment). At V g = 0, 
this peak is much less pronounced then at lower values 
of V g . This can also be seen in the dI/dV s d line traces 
of Fig.|3|>c and Fig.0J>c (solid curves). These observa- 
tions are similar to conventional S-N-S structures, such as 
atomic-sized break junctions [bH ]: for large transparen- 
cies of the junction a peak is observed at V s d — but no 
structure at 2A while for small transparencies a gap is ob- 
served around V s d = and a large peak at 2A marks the 
onset of quasi-particle tunneling. For these systems, the 
transparency depends on the atomic arrangement of the 
junction, here the effective transparency can be tuned by 
moving the level position of a single-electron state with 
a gate electrode. The effective transparency is large if 
the level is aligned with the Fermi levels of the leads (on 
resonance) and it is small otherwise (off resonance) . 

The subharmonic gap structure is clearly visible in the 
measured data of Fig. |3] and has a similar gate- voltage 
dependence as in the model calculation of Fig. 01 How- 
ever, there are several differences. The most dramatic 
one is the pronounced peak at (V g ,V s d) = (0, A/e) in 



the measurement (Fig. |3;). Because this position corre- 
sponds to an even MAR cycle it should be absent based 
on our previous consideration (see Fig. |2i) . 

Let us compare theory and experiment by focussing 
onto the dI/dV s d lines traces shown in Fig. ^p-c and 
Fig. 0J>c (solid curves). In panels b) the dot level is 
off resonance, while it is on resonance in panels c). For 
the former case, experiment and theory agree fairly well. 
The differences in MAR structure between model and 
experiment are much more pronounced at the resonance 
position. Whereas the experiment (Fig. 0;) reveals pro- 
nounced peaks at A/2, A and 2A, the calculated dI/dV s d 
(Fig. solid) reveals fine structure for small V s d and 
pronounced peaks at 2 A/3 and 2 A. According to our 
previous discussion dI/dV s d should indeed show a pro- 
nounced peak at V s d — 2A/3e, if the dot level is centered 
in the middle, i.e. for ip g = at point C in Fig. |2K 
and Fig. It rather appears in the experiment that, 
contrary to expectations, the subgap feature at 2 A/3 is 
missing, while the 'forbidden' at A is present. Such be- 
havior would only be expected for very asymmetric junc- 
tions having C s /Cd <C 1 [lj|, which is not the case in the 
present work. 

There are different imaginable scenarios that may ac- 
count for the observed A peak and the lack of fine struc- 
ture around V s d = mV in the data of Fig. 121;. Inelastic 
scattering processes inside the dot, for example, would 
broaden and obscure higher-order MAR features. Other 
possible reasons may be found in a broadened BCS spec- 
tral density (the superconductor consists of a bilayer of 
Au/Al or a suppression of higher-order MAR due 
to the on-site Coulomb repulsion. 

In a phenomenological approach, we may try to ac- 
count for the additional broadening by manually intro- 
ducing larger bare couplings T S) d- Many curves with 
varying parameters were calculated of which a represen- 
tative set is displayed in Fig. 0J)-c (dashed curves) cor- 
responding to relatively large dot-electrode couplings of 
r s = 2.5 A and Td = 3.5 A. For the off-resonance po- 
sition (Fig. 4b) the main effect of the larger T is the 
increased magnitude of dI/dV s d- In contrast, the MAR 
structures significantly changes for the resonance position 
(Fig. 4c). Remarkably, at large coupling T, peaks now ap- 
pear at 2 A, A and A/2. These peaks do not originate 
from the resonant level, but from the two neighboring 
ones which are off resonance (the dot levels are spaced 
by AE w 6.5 A). Though the agreement is now reason- 
able, there is one remaining problem. We were unable to 
reproduce the relative peak height between the 2A and 
A peaks. Using any reasonable set of parameters, the 
2A peaks is always larger than the A peak in the model, 
while it is the opposite in the experiment. We emphasize 
that the model does not take into account interaction 
and correlations. Since a Kondo resonance is observed 
in the normal state, which need not be suppressed in the 
superconducting state ^3] > this may be the origin of the 
discrepancy. The Kondo resonance changes the spectral 
density in the leads by adding spectral weight to the cen- 
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tcr of the gap and removing spectral weight from the gap 
edges. The former tends to enhance the A peak, while 
the latter tends to suppress the 2A one. This expla- 
nation is attractive, but more work both in theory and 
experiment is needed to substantiate it. 

In conclusion, we have investigated the non-linear con- 
ductance characteristics of a quantum dot coupled to su- 
perconducting electrodes. We find a strong dependence 
of the MAR structure on the level position of the single- 
electron states. The experimental data is compared with 
a theoretical model, assuming a BCS density-of-states in 
the electrodes and an interaction- free dot. Reasonable 
agreement is possible, if the tunneling coupling to the 
leads is enhanced by a factor ~ 2 in the model as com- 



pared to the experimental value. There are additional 
subtle differences which point to the importance of inter- 
action and exchange correlations. 
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